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Abstract 

We investigated the dissipative dynamics of quantum discord for correlated qubits 
under Markovian environments. The basic idea in the present scheme is that quantum 
discord is more general, and possibly more robust and fundamental, than entanglement. 
We provide three initially correlated qubits in pure Greenberger-Horne-Zeilinger (GHZ) 
or W state and analyse the time evolution of the quantum discord under various dissi- 
pative channels such as: Pauli channels a x , cr y , and a z , as well as depolarising channels. 
Surprisingly, we find that under the action of Pauli channel cr x , the quantum discord of 
GHZ state is not affected by decoherence. For the remaining dissipative channels, the 
W state is more robust than the GHZ state against decoherence. Moreover, we compare 
the dynamics of entanglement with that of the quantum discord under the conditions in 
which disentanglement occurs and show that quantum discord is more robust than en- 
tanglement except for phase flip coupling of the three qubits system to the environment. 
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1 Introduction 

One of the most remarkable properties of quantum mechanics is represented by the quantum 
correlation. Entanglement, which is a prominent feature of quantum correlation plays an im- 
portant role in quantum computing and informational processing [1, 2, 3, 4]. Recently, it has 
been perceived that entanglement is not the only kind of quantum correlation. In this con- 
text, others suitable measures of quantum correlations such as quantum discord [5], quantum 
deficit [6, 7, 8], quantumness of correlations [9] and quantum dissonance [10] have been pro- 
posed. Among them, quantum discord as a measure that is based on the difference between 
two quantum extensions of classically equivalent concepts has received considerable attention. 
This measure which quantifies the all nonclassical correlations between parts of a quantum sys- 
tem, actually supplements the measure of entanglement. For pure entangled states quantum 
discord coincides with the entropy of entanglement. It can also be nonzero for some mixed sep- 
arable state. It is worth mentioning that quantum discord is related to other concepts such as 
Maxwell's demons [11, 12], quantum phase transitions [13, 14], completely positive maps [15], 
and relative entropy [16]. Moreover, the characteristics of quantum discord have been studied 
in some physical models and information processing. It was shown that quantum discord can 
be considered as a more universal quantum resource than quantum entanglement in some sense. 
It offers new prospects for quantum information processing [17, 18, 19, 20, 21]. Studying of 
the quantum discord evolution exposed to noisy environments has led to the surprising result 
that it may obviously differ from entanglement evolution. In fact, the unavoidable interaction 
of the systems with their environment, based on completely positive quantum dynamical semi 
groups, can be modelled by means of the noisy quantum channels. Quantum channels are com- 
pletely positive and trace preserving maps between spaces of operators. The study of those 
can be broadly divided into Markovian and non-Markovian channels depend on the interaction 
of the system and environment. Markovian channels describes memoryless environments. The 
prototype of it is given by a quantum dynamical semi group, that is by solving a master equa- 
tion for the reduced density matrix with Lindblad structure [22, 23, 24]. For non-Markovian 
channels, environmental memory plays an important role, so the master equations describing 
their dynamics are often complicated integro-differential equations which are rarely exactly 
solvable [23, 24]. 

When a system of qubits with quantum correlation is exposed to noisy channels disentan- 
glement can occur suddenly, but the quantum discord mostly decays in the asymptotic time 
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[25, 26, 29, 27, 28]. This points to a controversial fact that quantum discord may be more 
robust against decoherence than entanglement. Hence quantum algorithms that are based only 
on the quantum discord correlations can be more robust than those based on entanglement. 
The aim of this paper is to illustrate the mentioned subject for a system of three-qubit which 
is initially prepared in pure state by Greenberger-Horne-Zeilinger (GHZ) [30] or W [31] state 
as 

\^ GHZ ) = ^(|000) + |in)) 

W) = |(|100) + |010> + v / 2|001». (1-1) 

Jung et al. [32] showed that these pure states will be mixed due to transmission through some 
of the common channels for qubits. Moreover, it was shown that under these noisy channels 
the sudden death of entanglement of three qubits occurs in a finite time [33]. The question is, 
what happens to quantum discord under the same conditions in which disentanglement can 
occur? 

In reply to this question we lead to remarkable results. The decoherence induced by the Pauli 
channel a X) can not affect the quantum discord of GHZ state contrary to the W state. We 
also show that the W state lose less quantum discord than the GHZ state due to transmission 
through the Pauli channels a y , a z and the depolarising channel. Comparison of entanglement 
and quantum discord demonstrates that for the Pauli channels a x and a y and the depolarising 
channel, quantum discord is more robust against decoherence than entanglement. 

The organisation of this paper is as follows: Section 2 is devoted to the necessary theoretical 
background to describe the time evolution of the system and introduces the global quantum 
discord. Evolution of quantum discord in transmission through the Pauli and the depolarising 
channels for the GHZ and W states is calculated in section 3 and 4, respectively. Finally, we 
summarise our results. 

2 Method 

2.1 Time evolution of states under Markovian channels 

In an open quantum system the prototype of a Markov process is given by a quantum dynamical 
semigroup of a completely positive and trace preserving map = exp [»££], t > 0. In this 
case, quantum evolution of the system is given by the solution of a Markovian master equation 
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with Lindblad structure for the reduced density matrix [22], 
j t p(t) = £p(t) 

= -z[H s ,p(t)] + J2^p( t ) L t-\{ L t^p( t )}- I 2 " 2 ) 

i 

For any Pauli channel a a [a = x,y,z), the decoherence dynamic is described by Lindblad 

A . 

operators La^o. = ,JKA~a. o"q 3 , j = 1,2,3, which act independently upon the j-th qubit. In 
these operators, ot 3 denote the Pauli matrices of the j-th qubit and the constants ka^o. are 
relaxation rates. For depolarising channel nine of these operators are needed. Here we assume 
that the Hamiltonian of the system is zero H s = 0, and the strength of the coupling between 
each of the qubits and channels is equal. Recently, Jung et al. [32] analysed time evolution of 
three qubit GHZ and W states in the presence of noisy channels. Here we use their results. 



2.2 Quantum discord 

For a bipartite system AB quantum discord is given by [5] 

D(p AB ) = I(p AB ) - C(p AB ), (2-3) 

where I(p AB ) = S(p A ) + S(p B ) — S(p AB ) ) is quantum mutual information which includes the 
total correlation between A and B. The last section on the right represents classical correlation 
C(p AB ) = max {Uk} [S(p A ) - S(p AB \{U k })] with S(p) = -Tr[p log 2 p] as the von-Neumann 
entropy. Notice that the maximum is taken over the set of projective measurements {n^} [35]. 
By definition the conditional density operator p AB = ^{(I A <8> nf )p AB (I A ® nf )} with p k = 
Tr[(I A ®U B )p AB ] as the probability of obtaining the outcome k. We can define the conditional 
entropy of A as S(p AB \{Il k }) = J2kPkS(p A ). This entropy includes the knowledge of subsystem 
B, with p A = Tr B [p AB ] and S(p A ) = S{p AB ). It has been shown that D(p AB ) > with the 
equal sign, only for classical correlation [36]. 

Very recently, Rulli et al. [34] have proposed a global measure of quantum discord based on a 
systematic extension of the bipartite quantum discord. Global quantum discord (GQD) which 
satisfy the basic requirements of a correlation function, for an arbitrary multipartite state 
pA x ...A N uric } er a se t Q f i oca i measurement {H^ 1 <g> ... <g> H^ N } is defined as 

N 

D(p A - A ») =™[5(/ 1 "^||$(/ 1 "^))-5:S(^||$ i (^))]. (2-4) 

Where ^j(p^) = E l Tlf j p A mf j and $(p^-^) = U k p A '- A ^U k with n fc = n^ 1 ® ...<8>n^ 
and k denoting the index string {ji-.-Jn)- We could eliminate dependence on measurement by 
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minimising the set of projectors {H^ 1 , ...,11^}. 

With these remarks about the global quantum discord (4), one can describe the time evolution 
of the quantum discord for three-qubit GHZ and W states when they are passed through a 
noisy channel. By selecting a set of von-Neumann measurements as 



Aj _( cos 2 (|) e^cos(|)sin(|) 



Ve-^cos(f)sin(f) sin^) )' 

jrA, _ ( sin 2 (|) cos(|) sin(|) \ 

2 " Ue^cos(|)sin(|) cos 2 (|) J' ^ 

the quantum discord for p = p GHZ (t) or p w (t) can be obtained from 

D(p) = min [S(p\Mp)) ~ E S(p^<M/0)]- (2-6) 
Where ^ G [0, 7r) and G [0, 27r), for j = 1,2,3, are azimuthal and polar angles, respectively. 



3 Quantum discord of three-qubit system with initial 
GHZ state under Markovian noise channel 



3.1 Pauli channel a x 

When a three-qubit system with initial GHZ state is coupled to a shift-flip noise channel, 
in which each qubit is coupled to its own channel, the time evolution is obtained by the 
solution of the master equation (2) with Lindbald operators Lj :X = yjK^ x Ox 3 , (j = 1, 2, 3). For 
this coupling, the master equation reduces to 8 diagonal and 28 off-diagonal coupled linear 
differential equations. So after transmission of the GHZ through the Pauli channel a x the 
density matrix is given by [32] 
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(3-7) 
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with 

a + = l + 3e' 4Kt , 

a_ = l-e~ 4Kt . (3-8) 
For this case, the lower bound to the concurrence is [33] 

r,{p G x HZ {t)) = e-^. (3-9) 

In order to find the analytical expression for the quantum discord with p GHZ (t) we must 
consider equation (6). By tracing out two qubits, the one qubit density matrices representing 
the individual subsystems are proportional to the identity operator 

pi M \t) = P i A Ht) = P^\t) = ^±^z 7 . (3.10) 
Therefore, we have S(p^\t)\\^j(p^\t))) — (j — 1,2,3) and the equation(6) reduces to 

D{ P G : HZ {t)) = ^n{s(p^ z (t)\Mp GHZ (t)m 

{0j,(pj} 

= ™n{S(<S>(p GHZ (t)))-S(p GHZ (t))}- (3-H) 

{0j,ipj} 

The entropy S(p GHZ (t)) can be obtained as 

S{p G x HZ {t)) = 2 - ^ log 2 q_ - ^ log 2 a + . (3-12) 

In order to obtain the maximum classical correlation among the parts of p GHZ (t), by varying 
the angles 9j and (pj, we must find the measurement bases that minimise D(p^ HZ (t)). After 
some calculation we have perceived that, by adopting local measurements in the a z eigenbases 
for each particle, the value of quantum discord will be minimised. So the von- Neumann entropy, 
after completing such measurements, can be written as 

S(Hp GHZ (t))) = 3 - ^ log 2 a_ - 2± l 0g2 a+ . (3-13) 

By substituting equations (12) and (13) into equation (11) quantum discord is readily found 
to be 

D{p G x HZ {t)) = 1. (3-14) 

As can be seen, the Pauli channel a x can not change the quantum discord of the three qubit 
systems with the initial GHZ state. We have plotted the dynamic evolution of the quantum 
discord for density matrix (7) versus the dimensionless scaled time nt in Fig. 1 (dashed violet 

)• 
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3.2 Pauli channel cr. 



If the three-qubit GHZ state are transmitted through the Pauli channel a y , its density matrix 
takes the following form [32] 



P° HZ (t) 
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(3-15) 



where a± are given in equation (8) and, ft and ft are defined as 



2ret , — 6ret 
1~ C , 



ft = 3e 

ft EE e -2«t_ e -6 K t 



(3-16) 



For this matrix, the entanglement vanishes after some finite time due to the condition [33] 

1 



?3 (/# (*)) = max{0, j(3e" 2Ki + e" 4Kt + e~ bKt - 1)}. 



(3-17) 



The reduced density matrices of subsystem of (15), is the same as equation (10). Hence, the 
quantum discord becomes 



(3-18) 



D(p^(t)) = min x {S(<S>(p^(t))) - S(p^(t))}- 
When the projective measurement in the eigenprojectors of a z is performed on any one of the 

rXj. 

V 



remaining qubits of p^ HZ (t), the minimum quantum discord is obtained as 



log 2 (a+ - ft) + 



log 2 («+ + ft) 



3(a_-ft) 3(a_ + ft), , a , 
+ g log 2 (a_-ft) + log 2 (a_ + ft) 

~ -j- log 2 (a+) - — log 2 (a_). 



(3-19) 



Quantum discord is reduced due to the Pauli noisy channel a y and disappears with low speed, 
as shown in Fig. 1 (solid blue ). 
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3.3 Pauli channel a z 

Transmission of the GHZ state through the Pauli channel a z result in the master equation (2) 
reduces to 8 diagonal and 28 off-diagonal first order differential equations with a simply trivial 
solution. So p GHZ (0) = \GHZ){GHZ\ is evolved into [32] 

p G z HZ {t) = i(|000)(000| + |111)(111|) + ie- 6K *(|000)(f f f | + |000)(111|). 

(3-20) 

For this mixed state, the lower bound to concurrence is a mono-exponential function of time 

r 3 (p GHZ (t)) = e~ 6Kt . (3-21) 

Let us now focus on the quantum discord for p GHZ {t) as defined by (6). After tracing out two 
qubits, the three reduce density matrices are equal, given by p z Al \t) = p z A ' 2 \t) = p[ As \t) = |. 
Hence, we have S , (pi A, ' ) (i)||$ J -(pi A, ' ) (t))) = (J = 1,2,3) and the equation (6) reduces to 

D( P GHZ (t)) = mm {Sm P GHZ (t))) - S(p GHZ (t))}- (3-22) 

The von-Neumann entropy of p GHZ (t) is 

S{ P G z HZ {t)) = {1 " 1± f^ log 2 (l + e~ 6 -) - log 2 (l - e" 6 -)}. 

(3-23) 

We have found that for the density matrix (20), measurements in the eigenprojectors of a J z , 
that is boundary values ifj = and 9j = 0, minimise &(p GHZ (t)) as 

*(p GHZ (t)) = \{\ + ++)(+ + +| + | - — )( 1}, (3-24) 

which cause S(&(p GHZ (i))) = I. In these circumstances, the quantum discord is expressed as 

D{p G z HZ {t)) = i±|— log 2 (l + e-^) + log 2 (l - e-^). (3-25) 

In Fig. 1, we have plotted the dynamic evolution of the quantum discord for the density matrix 
(20) as a function of nt (solid red). It can be seen that the quantum discord decreases from 
its maximal value D(p GHZ (t)) = 1 and vanishes after some finite time. 
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3.4 Depolarising channel 



For depolarising noise which is described by nine Lindblad operators, Lj iZ , Lj X and Lj y (j = 
1, 2, 3) the state of three-qubit system that were initially described by the GHZ state replaces 



P° d HZ (t) = ^ 
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(3-26) 



Since pp\t) 



a++3a_ 



/, then S(p^\t)\\^(p^>(t))) 



(3-27) 

0. The quantum discord, due to the 



condition of depolarizing channel may be written as 

D( P r Z (t)) = min{S($(pf^))) - S{p G d HZ {t))}- 



(3-28) 



We find that measurements in the eigenprojectors of a{ minimise the quantum discord. After 
some algebraic calculation it is found that 

(«+ + 7) 



D{ P r z (t)) 



8 



■log 2 (a+ + 7) 



+ (a+ g 7) log 2 («+ ~ 7) ~ ^ log 2 (a+)- 



(3-29) 



The time evolution of quantum discord (6) as a function of nt in the case of the depolarising 
channel is plotted in Fig. 1 (solid black). Due to the noisy channel, quantum discord decays 



from its maximum value D(pf HZ (0)) = 1 and disappears after some finite time. 



As we have pointed out previously, transmission of the GHZ state through the Pauli channel 
a x contrary to other channels, could not disturb the quantum discord. We also observed that 
the quantum discord of the initial state due to transmission of the GHZ state through the Pauli 
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channel a y decayed less quickly than a z and the depolarising channel, as shown in Fig. 1. In 
order to comparse quantum discord and entanglement we have used the results of reference 
[33]. In Fig. 2, the time dependent entanglement of the GHZ state under the noisy channels is 
shown. In all cases, one can see decrease of the entanglement due to the noise of the channels. 
Clearly, under the dissipative dynamics considered here, except for phase flip, coupling of the 
three qubits system to the environment, quantum discord is more robust than entanglement. 
In the next section we will discuss the effects of noisy channels when we prepare three qubits 
with initial W state. 



4 Quantum discord of three-qubit systems with initial 
W state under Markovian noise channel 



4.1 Pauli channels a x and a y 

If the three-qubit system is initially prepared in the W state, a similar analysis as in the 
previous section to compute the quantum discord can be made. Under the bit flip or bit- 
phase flip coupling to environment the density matrix of three-qubit W state at time t has the 
following analytical expression 



(*) 





/ 2a 2 








±V2a 2 





±V2a 2 


±a 2 


\ 







2ai 















±a 3 







>/2<*l 


2/3+ 





oil 








±V2a 3 


1 


±V2a 2 








2/3_ 







v^2a4 





16 












2/3+ 








±V2a 3 




±y/2a 2 








CH4 





2/3- 









±«2 








V^2«4 





V / 2a4 


2a 4 







V o 


±a 3 


±V2a 3 





±V2a 3 








2a 3 J 



(4-30) 

where the + sign refers to the a x and — to the a y channel, respectively. The density matrix 
elements are given by 



1 + e~ 2Kt + e~ iKt + e" 6K * 



a 2 = l + e- 2Kt 



-Ant 



—Gut 



a 3 = 1 

«4 = 1 
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0± = l±e- bKt . (4-31) 

It worth mentioning that since two density matrices p±(t) have the same structure of matrix 
elements, we expect that the quantum discord dynamic of these density matrices coincide. 
After tracing out any two qubits, the reduced density matrices are given by 

(a ) 1 /2 + e~ 2K * \ 

To find the measurement bases that minimise quantum discord, we used the same procedure 
as in the previous section. We perceived that the maximum classical correlation for any 
reduced density matrix (32) was obtained by completing the projective measurements in the 
eigenprojectors of a z . It leads to S(p { ± j \t) \\<f>j(p± j \t))) = 0. The minimum of D(p™(t)) 
and D(pY(t)) is obtained after the two project measurements in the eigneprojectors of a x and 
a z on p^{t) and p^(t), respectively. We do not show the achieved analytic expressions for 
the quantum discord here because its form are not compact. We have plotted the quantum 
discord for density matrix (30) as a function of the dimensionless scaled time nt at Fig. 3 
(dashed violet). Here, we can see that quantum discord due to the transmission through a x 
and tjy channels, as expected, coincide and decrease at a low rate of speed from its initial value 
D(p±(0)) = 1.5 and the asymptotic limit select the exact value 0.813. 



4.2 Pauli channel o z 

After transmission through the Pauli channel a Z) the time evolution of the W state is described 
by [32] 
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(4-33) 



The density matrices representing each of the one qubit subsystem of (33) are given by 
p^ 3 \t) = |{3|0)(0| + |1)(1|} (j = 1,2,3). By using the same procedure as above to find 
the measurement bases that minimise D(pY(t)), we find that by performing the projective 
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measurement in the eigenprojectors of a J z quantum discord reaches its lowest point. Thus one 
can verify that S(p^ 3 \t)\\<&j(pz^ 3 \t))) = and S($(pY(t))) = §• The von-Neumann entropy 
of pf if) density matrix is found as 

S(pY(t)) = 1(11 + e -^)- 1(1 - e -^)tog 2 (l-e-^) 

- -{(3 + e~ iKt - Vl - 2e~ 4K * + 17e- 8K ') 
log 2 (3 + e~ 4K * - VI - 2e" 4K * + 17e~ 8Kt ) 

+ (3 + e~ 4Kt + Vl - 2e' iKt + 17e~ 8Kt ) 

log 2 (3 + e" 4Ki + Vl - 2e- 4K * + 17e- 8K ')}- (4-34) 

Therefore, the expression for quantum discord D(pY{t)) is given by 

D( P Y(t)) = ^(5 + e-^) + l(l-e- 4Kt )log 2 (l-e- 4Kt ) 

+ i{(3 + e- 4Kt - Vl - 2e" 4Kt + 17e- 8K *) 
log 2 (3 + e" 4K * - Vl - 2e" 4Ki + 17e- 8Kt ) 

+ (3 + e" 4Ki + Vl - 2e- 4Kt + 17e" 8 «*) 

log 2 (3 + e~ 4Kt + Vl-2e- 4K * + 17e- 8K ')} . (4-35) 

With the above channels, quantum discord of the selected system with the initial W state 
disappear under Pauli channel a z after a finite time. This is displayed in Fig. 3 (solid red). 

4.3 Depolarising channel 

When three-qubite W state is transmitted through the depolarising channel, its time evolution 
is described by [32] 
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(4-36) 
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(4-37) 

For density matrix (36), the single qubit density matrix is given by 

(A) 1 /2 + e" 4K * \ 
^(*) = 4( o 2-e--j J = W - (4 " 38) 

Our results show that the minimum quantum discord for the depolarising coupling of the three- 
qubit W state is attained at boundary values <fj = and 9j = 0, that is for a local measurement 
along the eigenstates of Pauli matrix a z . The state of the single qubit is not induced by such 
measurement, therefore S(p { d Aj) (t^^p^ (t))) = 0. The von-Neumann entropy of pY (t) after 
measurement is given explicitly by 

S($(p%(t))) = -{18-ailog 2 ai -a 2 log 2 a 2 

- a 3 log 2 a 3 - a 4 log 2 a 4 - P+ log 2 f)+ - /L log 2 /L } . (4-39) 

In this case, the time evolution of quantum discord is obtained as 

D( P Y(t)) = i(3-e- 8K ')-i{^log 2 ^ + + ^_log 2 ^_ 

+ oli log 2 oli + 2a 2 log 2 a 2 + 2a 3 log 2 a 3 + a A log 2 fi 4 } 

+ \{{fi- ~ 7-) Iog 2 ()9- - 7-) + 0+ ~ 7+) log 2 (^+ " 7+)} 

+ ^{(/3 + + 7+ + «i - + 7+(2^+ + 17%) - 2d^0~ + + 7+) + a?) 

log 2 (/3+ + 7+ + «i - V/3+ + 7+(2/3 + + 177+) ~ 2<*i0+ + 7+) + &i) 
+ 0+ + 7+ + ai + + 7+ (2/3+ + 177+) ~ 2ai(/?+ + 7+) + jg) 

log 2 (/3+ + 7+ + «i WW± 7+(2/3+ + 177+) - 2ai(/3+ + 7+) + 
+ 0_ + 7- + a 4 - V/fe + 7- (2/9- + 177_) - 2a 4 (^_ + 7-) + g) 



log 2 (/3_ + 7_ + S 4 - V/^ 2 + 7- (2)9- + 177-) - 2<5 4 (/3_ + 7-) + «1) 



+ (/?_ + 7_ + a 4 + V/? 2 + 7- (2/3- + 177-) - 2a 4 (/?_ + 7 _) + a 2 ) 



log 2 (/3_ + 7_ + a 4 + V/3 2 + 7-(2/3_ + 17%) - 2%(% + 7.) +a\)}. 



(4-40) 



The responses of quantum discord with the initial W state to the different noises as functions 
of Kt are listed at Fig. 3. As in the case of the GHZ state, the quantum discord for the W 
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state decays due to the noisy channels. When a three-qubit system with initial W or GHZ 
state is coupled to a depolarising noise channel decay of the quantum discord occurs so fast. 
This means that the depolarising channel has the most destructive influence on the quantum 
discord. However, the quantum discord decreases due to the transmission through the Pauli 
channels a x or a y , but it catches the exact value 0.813 in the long term as shown in Fig. 3. 
This figure also exhibits the vanishing of the quantum discord for the initial W state under the 
phase flip channel. These results show that GHZ-type quantum discord is more fragile under 
certain types of environment coupling which can be modelled by means of the Pauli channels 
o y and o z as well as the depolarising channel as compared to the quantum discord of the W 
state. However, in the Pauli channel a x the situation becomes completely reversed. 
In order to show the difference between the entanglement and quantum discord under the 
various dissipative channels, by using the results of [33] we have plotted the evolution of 
entanglement of the W state in Fig. 4. As in the case of the GHZ state, the entanglement 
for the W state decays exponentially due to the noisy channels. Observe that under the 
dissipative dynamics considered, except for the phase flip coupling of the three qubits system 
to the environment, discord is more robust than entanglement. 

5 SUMMARY AND CONCLUSION 

To sum up, we prepared the three-qubit system with the initial state formed by GHZ and W 
states. We then studied the dynamics of quantum discord under interaction with indepen- 
dent Markovian environments which can be modelled by means of the various noisy channels, 
namely, the Pauli channels a x , o y and o z as well as the depolarising channel. It was clarified 
that coupling with the bit flip channel could not disorder the quantum discord of three-qubit 
GHZ state unlike in the W state. In other words, for this noisy channel, GHZ-type quantum 
discord is always more robust than quantum discord of the W state. In the Pauli channels 
dy, o z and depolarising channel, the W state preserved more quantum discord than the GHZ 
state in the long term. 

Also, we observed that under the dissipative Markovian dynamics considered, except for the 
Pauli channel a z quantum discord is more robust than entanglement. This points to the fact 
that quantum discord is another kind of quantum correlation different from entanglement and 
the absence of entanglement does not necessarily indicate the absence of quantum correlations. 
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Figure Caption 

• FIG. I. (Colour online) The quantum discord (6) for the three-qubit system with the 
initial GHZ state as a function of nt, if transmitted through Markovian channels: Pauli 
channel o x (dashed violet), o y (solid blue), a z (solid red) and depolarising channel (solid 
black) . 

• FIG. 2. (Colour online) The entanglement for the three-qubit system with the initial 
GHZ state as a function of Kt, if transmitted through Markovian channels: Pauli channel 
a x (dashed violet), a y (solid blue), a z (solid red) and depolarising channel (solid black). 

• FIG. 3. (Colour online) The quantum discord (6) for the three-qubit system with the 
initial W state as a function of nt, if transmitted through Markovian channels: Pauli 
channel o x and a y (dashed violet), a z (solid red) and depolarising channel (solid black). 

• FIG. 4. (Colour online) The entanglement for the three-qubit system with the initial W 
state as a function of nt, if transmitted through Markovian channels: Pauli channel a x 
(dashed violet), a y (solid blue), a z (solid red) and depolarising channel (solid black). 



